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We apply strong-coupling expansion techniques to finite-temperature 
lattice pure gauge theory, obtaining dimensionally reduced Zjv-symmetric 
effective theories. The analytic mappings between the effective couplings 
and the original one, viz. j3, allows to estimate the transition point f3c 
of the 4D theory for a large range of the imaginary-time extent Nr of 
the lattice. We study the models for SU{3) via Monte Carlo simulation, 
finding satisfactory agreement with the critical point of the original theories 
especially at low Nr- 

PACS numbers: 11.15.Ha, 12.38.Aw 

1. Introduction and theoretical setting 

Among the possible approaches to the study of Lattice QCD, effective 
theories play an important role, sometimes opening the way to subjects 
otherwise inaccessible, giving a deeper understanding of the physics at play, 
or even simply reducing the computational efforts involved. A desirable 
condition is that the effective theory is motivated by first principles and 
retains the original symmetries. The history of QCD effective theories is 
rather long; this work aims at providing the final results for the strong- 
coupling approach presented in 

The effective theories considered in this work re-express the partition 
function of a (3 -|- l)-dimensional Yang-Mills SU{N) lattice theory (all ex- 
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plicit calculations refer to = 3) as a three-dimensional model with com- 
plex numbers as per-site degrees of freedom, the (traced) Polyakov loops in 
the 4D system = TrJlfli Uo{x,t): 

Z = l[dU]exp[-^Y.^eTrUa) ^ J [dL,]e^'''^^^'^^+^^^^W+-) . (i) 

The effective models will exhibit various spin-like interaction terms Sn, 
each with a specific coupling A„ which is a function of the imaginary-time ex- 
tent Nr of the 4D lattice and its (bare) coupling /3. The strong-coupling se- 
ries for the A„(/3, Nr) employs a character expansion and then the moment- 
cumulant formalism [1112]; here we only report the final formulae. 

Of the (infinitely many) interaction terms, we consider the three featur- 
ing the lowest order in /3 (or in n = a/(/3) ~ /3/18 + 0(/3^), see e. g. [2]): 
nearest- and next-to- nearest-neighbour fundamental representation (Ai5i, 
and nearest-neighbour adjoint interaction (XaSa)- We study three 
different effective theories: one with only the Ai interaction, one with 
(Ai, A2), and the last with the (Ai, A^) terms. We parametrise L as0 

LAOx,(t>x) = e'' + e'^ + e-'^'+^^ , -tt < 9,,^, < +7r ; (2) 

f dW, ^ [ "d(/.,(27-18|L,|2 + 8ReL3 - |L,|^) .(3) 

JSUi'i) J-T7 J-TT V ' 

=exp(y,) 

The effective theories studied numerically are given by0 Zn) = Zn 



(1) = ^^(1,2)^2=0, 



%,2) = n/d^- / d^xll'^Myx) n (l + 2AiReL,L*)n(l + 2A2ReLfcLn , 

^ ^ <i,j> lk,l] 

= n / / d^xll'^Myx) n (1 + 2AiReL,L*) • 

^ X <i,j> 

■ n [l + Xa{\Lm\'-l)i\Ln\'-l)] , (4) 

<m,n> 

In [T] we presented the mappings A„ -H- /3 for even A^,-; here we provide in 
addition A'',- = 1, 3: Xi{u; 1) = u, A2(n; 1) = and Xa{u; 1) = v; for Nr = 3 
the Xa map follows the general formula, while: 

287 

Ai(n; 3) = exp[3(4u'' + 12^^ - Uu^ - 36n'^ + —u^ + 

1851 q 932917 m,, 

H H u ) ; 

10 5120 ' 

A2(u; 3) = u^{Qu^ + 18u^ + 117u^) . (5) 



^ We note that the measure used in [1] contains an error corrected here. 
^ <i,j> denotes nearest-neighbours and [k,l] next-to-nearest-neighbours. 
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Fig. 1. One-coupling model for various system sizes. Left: behaviour of \L\ as 
a function of Ai and (inset) scatter plot of for a small and a large coupling. 
Middle: Binder cumulant Right: Susceptibility X\l\- The vertical line marks 

the phase transition. 



2. Numerical simulations and phase structure 

All three models were simulated on cubic systems with A*"^ sites and peri- 
odic boundary conditions via a Metropolis accept/reject procedure. Looking 
at the expressions in Eq. 21 one realises that at sufficiently high couplings a 
"sign problem" might occur for negative values of Re(LjL^): to take care of 
this, the simulations uses weights |(1 + 2AiReLjL|)| (and similarly for the 
other terms) and folds the sign into the observable, which is subsequently 
reweighted to get the correct answer. It turns out that, within the range 
of couplings of interest, this problem is very mild and the average sign of a 
configuration never drops below 0.999. 

One-coupling model. A first inspection of the distribution of at dif- 
ferent couplings confirms the existence of a phase transition at some finite 
a[^^ (fig- 1])- More quantitatively, we use as basic observable 

We measure the Binder cumulant and the susceptibility: 

%i = i-^|^ ; x\L\ = m')-m? ; (7) 

the minimum of the former and the maximum of the latter (see fig. [T]) 
are used as size-dependent criticality estimators Xb,x{^s), and a finite-size 
scaling fit is then attempted on both, with first-order scaling law 

X{Ns) = X^ll + bN^^ . (8) 
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Fig. 2. Scaling analysis for the one-coupling model. Left: data and fit to Eq. [5]for 
the two pseudo-criticality estimators. Middle: data and fit to Eq. [5] Right: data 
and fit to Eq. [SI with both estimates for the asymptotic values shown as horizontal 
lines. 



Data from systems of side Ng = 6, 8, 10, 12, 14 allowed for two indepen- 
dent and consistent estimates of the infinite-volume critical point, averaged 
to X^^l = 0.187885(30). Further evidence in support of the first-order nature 
of the transition comes from the scaling of the y-coordinates of the extrema 
of the observables in Eq. [3 we found for both (fig. [2D 

y^MNs) = y^.B + (const.) X N;^ , (9) 

with > and = 0.66277(7): the latter, as required, is lower than 2/3, 
and consistent with the estimate Ub oo ~ 0.6617(15) coming from locating 
the two maxima |-L|i,2 in the distribution of |L| at criticality [3]. 

Two-coupling model (Ai, A2). In the two-parameter cases, the (3c of the 
original gauge theory is found by intersecting the critical line of the models 
per se and one-dimensional manifolds encoding the original theory at A',-. 
Ten values of A2 < 0.01 were fixed and the approach of the previous case 
was applied to all of them: a polynomial fit to Ai = + aiX2 + a2A2 
gives the curve in fig. [3l compatible with the value for A2 = 0, with 
ai = {0.18787(2), -3.375(8), 12.8(7)}. 

Two-coupling model (Ai, Aq). The same procedure led to a parametri- 
sation of the critical line in the (Ai,Aa) plane, with 13 sampled values of 
Aq < 0.12. In this case the fitted function was Ai = cq + ciAq -|- C2A^ + csAj^ 
(fig. [3]), with coefficients a = {0.18783(8), -0.50(2), -6.9(7), 30.4(5.7)}. 

3. Results and conclusions 

The 4D data were taken from [4, 5], except the A',- = 1 critical point 
which was determined with a dedicated set of standard SU (3) Wilson action 
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Fig. 3. Two-coupling phase spaces: the critical line parametrisation is shown along 
with the measured datapoints and the iVT--specific manifolds, for the (1,2) and 
(l,a) models (resp. left and right). 
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1 


2.7030(040) 


2.78283(38) 




2.52906(613) 


2 


5.1000(500) 


5.18391(21) 


5.01735(36) 


5.00295(513) 


3 


5.5500(100) 


5.84878(11) 


5.73325(27) 


5.78014(181) 


4 


5.6925(002) 


6.09871(07) 


6.05229(11) 


6.07479(056) 


6 


5.8941(005) 


6.32625(04) 


6.32399(03) 


6.32250(011) 


8 


6.0010(250) 


6.43045(03) 


6.43033(02) 


6.42971(007) 


10 


6.1600(070) 


6.49010(02) 


6.49008(02) 


6.48991(006) 


12 


6.2680(120) 


6.52875(02) 


6.52874(01) 


6.52869(005) 


14 


6.3830(100) 


6.55584(02) 


6.55583(01) 


6.55580(004) 


16 


6.4500(500) 


6.57588(01) 


6.57587(01) 


6.57585(003) 



Table 1 . Critical couplings for various from different effective theories compared 
to the 4D Monte Carlo results. 



simulations. The values for l3c{N-j-) are summarised in Tabled] and plotted 
in fig. m along with a plot of the ratio between the effective-theory result 
and the full 4D Yang-Mills outcome. 

The discrepancy never exceeds ~ 7-8%; it is noteworthy that the one- 
coupling action seems to provide the best estimates at A'^^ = 1,2. At low 
Nt, where j3 is smaller, we expect the series expansion to show a better 
convergence; on the other hand, the particular shape of the mappings Ai(/3) 
is such that, as A'^,- increases, only the Ai coupling is important, which is 
the reason why the various results tend to converge one onto another. 

It must be stressed that the present results, able to reproduce the critical 
points with some accuracy, are obtained with a rather small computational 
effort (3D instead of 4D and complex numbers instead of matrices), mea- 
surable in the range of a few days with an ordinary desktop PC. This is in 
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Fig. 4. Comparison between the determinations of f3c{Nr) from the effective models 
and the 4D Monte Carlo results. Left: plot of Right: the ratio 



contrast with the "inverse Monte Carlo" approach 0, which requires sim- 
ulating the full theory in order to fix the coefficients; on the other hand, 
the latter technique is able to reproduce the theory on both sides of the 
transition. 

An extension of the present work is currently being carried on, with the 
introduction of massive fermions in the effective formulation via a hopping- 
parameter expansion. 
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